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ABSTRACT 

We present Spitzer/IRAC 4.5 transit photometry of GJ 3470 b, a Neptune-size planet orbiting a 
Ml. 5 dwarf star with a 3.3-day period recently discovered in the course of the HARPS M-dwarf survey. 
We refine the stellar parameters by employing purely empirical mass-luminosity and surface brightness 
relations constrained by our updated value for the mean stellar density, and additional information 
from new near-infrared spectroscopic observations. We derive a stellar mass of M^, = 0.539lg'o43 Mq 
and a radius of = 0.568toml i?0. We determine the host star of G J 3470 b to be metal-rich, 
with a metallicity of [Fe/H] = 4-0.20 ± 0.10 and an effective temperature of T^s = 3600 ± 100 K. The 
revised stellar parameters yield a planetary radius Rp = 4.83lg'2i -^9 that is 13% larger than the value 
previously reported in the literature. We find a planetary mass Mp = 13.9l};4Me that translates to 

a very low planetary density, pp = 0.72^° gcm^"^, which is 33% smaller than the original value. 
With a mean density half of that of GJ 436 b, GJ 3470 b is an example of a very low-density low-mass 
planet, similar to Kepler-lld, Kepler-lle, and Kepler-18c but orbiting a much brighter nearby star 
that is more conducive to follow-up studies. 

Subject headings: planetary systems - stars: individual (GJ3470) - techniques: photometric, spectro- 
scopic 



1. INTRODUCTION 

In the regime of low mass exoplancts only a handful 
of those known to periodically pass in front of their host 
stars have transits that are deep enough and orbit par- 
ent stars that are bright enough to make them amenable 
to extensive follow-up observations. The Kepler mission 
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has recently announced a harvest of more than 2,700 
planetary candidates identified since the launch of the 
spacecraft in 2009 (Batalha et al. 2012). About ten per- 
cent are Jupiter-size planets with radii between 0.7 and 
2.0 Jupiter radii, while more than 55% are Neptune- 
size planets with radii between 2 and 6 Earth radii. 
On the other hand, among the 241 confirmed transit- 
ing exoplanets (coming mainly from ground-based sur- 
veys), 62% are Jupiter-size planets with radii between 
0.7 and 2.0 Jupiter radii^^. It is now clear from Ke- 
pler and other studies that short-period Jupiter-size ob- 
jects make up a relatively small fraction of the exoplanet 
population (e.g., Howard et al. 2010; Wittenmyer et al. 
2011; Howard et al. 2012) . This stark contrast between 
confirmed exoplanets and the large underlying popula- 
tion glimpsed by Kepler has motivated intense efforts to- 
wards the characterization of smaller planets, in order to 
reach a comparable state of knowledge to what has been 
learned about the hot- Jupiter population. These efforts 
began already several years ago with the launch of a num- 
ber of ground-based projects dedicated to M-dwarf mon- 
itoring using both spectroscopy (e.g., the HARPS pro- 
gram; Bonfils et al. 2013) and photometry (e.g., MEarth; 
Nutzman & Charbonneau 2008). Planets orbiting M- 
dwarf stars offer the possibility to probe smaller plan- 
ets for a given transit depth, because of the favorable 
star-to-planet radius ratio. G J 436 b (Butler et al. 2004; 
Gillon etal. 2007) and G J 1214 b (Charbonneau et al. 
2009) are the smallest planets orbiting M stars with K 
magnitude brighter than nine, enabling detailed follow- 
up studies both from the ground and from space (e.g.. 

Source: http://www.exoplaiiets.org 
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Stevenson et al. 2010; Bean ct al. 2010). 

In the Neptune-mass range, thanks to its relatively 
large transit depth and host star brightness, GJ436b 
remains a "Rosetta stone" for our understanding of a 
whole class of exoplanets, shown to be ubiquitous in our 
Galaxy. With a mass 22 times that of the Earth and a 
radius 4 times larger than our home planet, G J 436 b 
has a relatively high density (pp = 1.69tg;^2 gcm~^; 
Torres et al. 2008), suggesting the presence of a mas- 
sive core made of silicates and/or ices. However, a 
H/He envelope is needed to reproduce its observed ra- 
dius (e.g., Figueira et al. 2009; Rogers & Seager 2010a). 
The improvement in the planetary radius of this ob- 
ject brought about by Spitzer observations placed signif- 
icant constraints on the range of possible compositions 
of GJ 436 b's interior. A key question that still needs to 
be addressed, however, is the extent to which GJ436b is 
representative of the entire exo-Neptune population. 

G J 3470 b is a new transiting Neptune-size planet dis- 
covered in the past year (Bonfils et al. 2012). It orbits 
a Kg = 7.99 mag, Ml. 5 dwarf with a period of 3.337 
days. With a published mass of 14.0 ± 1.7 Earth masses 
and a radius of 4.2 ±0.6 Earth radii (Bonfils et al. 2012), 
GJ 3470 b has a mean density pp = 1.07 ± 0.43 gcm~'^ 
that is significantly smaller than that of G J 436 b. The 
Kepler mission confirmed several of these so-called "low- 
density Neptunes". The first two were Kepler-lld and e 
(Lissauer et al. 2011), both belonging to the most pop- 
ulated transiting planet system known to date, and the 
third was Kepler-18c (Cochran et al. 2011), also a mem- 
ber of a multi-planet system. These objects represent 
the tip of the iceberg, as several hundred Neptune-size 
planet candidates have already been detected by Kepler 
and await confirmation. Unfortunately, most of these 
Kepler planets orbit faint stars and exhibit shallow tran- 
sit depths that render follow-up studies very challeng- 
ing, if not impractical altogether. Aside from Kepler, 
the ground-based survey HAT discovered the low-density 
Neptune HAT-P-26b (Hartman et al. 2011) which, until 
the discovery of G J 3470 b, represented the most promis- 
ing target for follow-up studies. However, as compared 
to G J 3470 b, the smaller planet-to-star area ratio cou- 
pled with the lower brightness of its larger Kl host star 
{K = 9.6) makes HAT-P-26b a less favorable target for 
follow-up studies. G J 3470 b therefore presents an ideal 
opportunity to investigate the internal structure, atmo- 
spheric composition, and possible formation pathways 
of low-density Neptune-size planets (e.g., Rogers et al. 
2011). 

All transit photometry available so far for GJ 3470 has 
been collected from the ground. While these time-series 
confirm the transiting nature of G J 3470 b, they do not 
precisely constrain the transit parameters, resulting in 
poorly determined planetary properties. We present in 
this paper the analysis of two transits of G J 3470 b ob- 
tained with the Spitzer Space Telescope at 4.5 fim in the 
frame of our DDT program (Demory et al. 2012b), sub- 
mitted shortly after GJ 3470 b's discovery. These data 
yield a significant refinement of the system parameters. 
The paper is organized as follows. Section 2 describes the 
observations and data reduction, while Section 3 presents 
the photometric and spectroscopic data analyses. Sec- 
tion 4 is dedicated to the stellar characterization, and 



the resulting planetary parameters are reported in Sec- 
tion 5. We discuss GJ 3470 b's internal structure and 
composition in Section 6. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Spitzer IRAC 4-5-pm photometry 

We observed two consecutive transits of G J 3470 b at 
4.5 /im using Spitzer^s InfraRed Array Camera (IRAC; 
Fazio et al. 2004). Observations took place on 11 June 
2012 and 15 June 2012 UTC as part of our DDT pro- 
gram PID 80261. For each transit we obtained 780 sets 
of 64 subarray frames each, with an exposure time of 
0.40 s per frame. Each Astronomical Observation Re- 
quest (AOR) lasted 6.5 hours, including 30 minutes over- 
head for the Pointing Calibration and Reference Sensor 
peak-up sequence. This step allowed GJ 3470 to be pre- 
cisely slewed on the position of maximum sensitivity on 
the 4.5/im channel subarray field-of-view (Ingalls et al. 
2012; Demory et al. 2012a). All data were processed by 
the Spitzer pipeline version S19.1.0, which produced the 
basic calibrated data necessary for our reductions. We 
first convert fluxes from the Spitzer units of specific in- 
tensity (MJy/sr) to photon counts, and transform the 
data timestamps from BJDutc to BJDtdb following 
Eastman et al. (2010). We then perform aperture pho- 
tometry on each subarray image using the APER routine 
from the IDL Astronomy User's Library We compute 
the stellar fluxes in aperture radii ranging between 1.8 
and 4.0 pixels, the best results being obtained with an 
aperture radius of 3 pixels. We use background annuli 
extending from 11 to 15.5 pixels from the Point Response 
Function center. For each block of 64 subarray images, 
we discard the discrepant values for the measurements 
of flux, background, and x-y centroid positions using a 
lO-cr median clipping for the four parameters. We then 
average the resulting values, the photometric errors being 
taken as the uncertainties on the average flux measure- 
ments. At this stage, a 50-ct clipping moving average 
is used on the resulting light curve to discard obviously 
discrepant subarray-averaged fluxes. Close examination 
of the resulting time-series reveals a sharp increase of 
the background and stellar fluxes, corresponding to the 
well-known "ramp" effect seen in other warm Spitzer ob- 
servations (see, e.g., Knutson et al. 2012). The raw pho- 
tometry for both AORs is shown on Figure 1. 

2.2. WIYN Speckle Observations 

We supplemented our G J 3470 b Spitzer photometry 
with speckle observations to explore the possibility of 
blended companions at close angular separations from 
GJ 3470. Speckle observations of G J 3470 were obtained 
at the WIYN 3.5-m telescope on 2 December 2012. The 
camera used was the Differential Speckle Survey Instru- 
ment, which is described by Horch et al. (2009). It is a 
dual-channel instrument that records images in two col- 
ors simultaneously. In the case of this observation, the 
filters used had center wavelengths of 692 nm and 880 nm, 
with filter widths of 40 and 50 nm, respectively. A speckle 
sequence of 3000 50-ms frames was taken on the target, 
followed by 1000 frames taken on a bright point source 
(HR3163) located near in the sky to GJ3470. These 

http: //idlastro . gsfc .nasa.gov/contents .html 
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Fig. 1. — Spitzer/IRAC 4.5/im photometry. Raw photometry 
from the two AORs is displayed after normalization. The optimal 
baseline model (logarithmic ramp model added to a second-order 
polynomial fit for the centroid position and a time-dependent trend 
(see Section 3) is superimposed for each AOR. in red. 

latter data are used as an estimate of the speckle trans- 
fer function for deconvolution in the reduction process. 
Reconstructed images are formed from the speckle data 
using the technique of bispectral analysis, which is de- 
scribed, e.g., by Horch et al. (2012). We then analyze 
the final images to determine the detection limits of faint 
companions near GJ 3470 using the technique described 
in the same paper. 

Figure 2 shows these detection limits based on the fi- 
nal diffraction-limited images in each filter. It is clear 
that there is no companion to the limit of our detection 
capabilities at a separation greater than 0'.'2. At 0'.'2 the 
limiting Am for the 692 nm image is 3.87 mag, and for the 
880 nm image it is 3.39 mag. Inside of this limit, as one 
approaches the central star, the limiting Am becomes 
smaller as the peaks and valleys of the reconstructed im- 
age get larger. In studying the two images, we find that 
none of the peaks near the central star are in the same 
position in both images, which is a good indication that 
they are probably not real stars but noise peaks. One of 
the advantages of the two independent channels in the 
instrument is to see if the positions of faint peaks match. 
We conclude that, to the limit of our detection at WIYN, 
there is no resolvable companion. 

2.3. Magellan/FIRE Near-infrared Spectroscopy 

Spectroscopic properties such as the effective temper- 
ature, Tcff, and metallicity, [Fe/H], needed to establish 
the physical parameters of the parent stars of transiting 
planets have usually been difficult to determine for M 
dwarfs. Several studies in the past year have presented 
calibrations of [Fe/H] or Toff in terms of easily measurable 
spectroscopic indices in the i/-band and ^iT-band regions 
that represent a significant advancement in the field. To- 
ward this goal, we obtained a near-infrared spectrum of 
GJ3470 on 12 November 2012 with the Folded-port In- 
frared Echellette (FIRE) spectrograph at the 6.5-m Mag- 
ellan Baade telescope. FIRE delivers R = 6000 spectra 
from 0.83 to 2.5 /zm in a single-object, cross-dispersed 
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Fig. 2. — Detection limit analysis of GJ 3470 based on speckle 
data. The squares and diamonds indicate the magnitude differ- 
ences of local maxima in each reconstructed image as a function 
of separation from the central star, and points are local minima. 
The curves indicate the b-a detection limit based on the statistics 
of the these maxima and minima. The dashed curve is the result 
for the 692 nm reconstructed image, and the dot-dash curve is the 
result for the 880 nm image. These curves indicate a greater than 
4-magnitude sensitivity in the limiting Am at most separations. 

setup (Simcoe et al. 2008). We used an exposure time 
of 8.5 minutes on GJ 347() and 105 seconds on the AOV 
telluric standard HD 58296. We reduce the data using 
fire's pipeline FIREHOSE, which employs the methods 
of Vaeca et al. (2003) for telluric correction. An internal 
ThAr lamp provides wavelength calibration of both the 
GJ 3470 and HD 58296 spectra. The signal-to-noise ra- 
tio in the reduced spectrum is > 200 in the K band, 
where H2O features near 2.21 and 2.26 /xm fall on order 
20 of fire's 21 cross-dispersed grating orders. Strong 
OH emission lines from the sky (which can introduce 
shot noise and residuals from sky subtraction) do not 
fall directly on these features. The measurement of var- 
ious spectroscopic indices from this FIRE observation is 
described in Sect. 3.2.1. 

3. DATA ANALYSIS 

3.1. Spitzer Photometry 

3.1.1. Baseline Model Selection 

We first perform an individual analysis of each Spitzer 
AOR to determine the optimal baseline model, which 
accounts for time- and position-dependent systematic ef- 
fects relevant to our IRAC 4.5 /im observations. We 
employ for this purpose our adaptive Markov-Chain 
Monte Carlo (MCMC) implementation described by 
Gillon et al. (2010). We test six baseline models of 
increasing complexity, and compare their Bayesian in- 
formation criteria (BIC; see, e.g., Gelman et al. 2003) 
to choose the baseline model that yields the highest 
marginal likelihood. We correct for the well-known 
"pixel-phase" effect using a second- to fourth-order x-y 
position-dependent polynomial, while the "ramp" is cor- 
rected using a second-order logarithmic model. We also 
check for time-dependent trends of instrumental and/or 
stellar origin by adding linear or quadratic functions of 
time to our baseline models. We additionally explore the 
correlation of the stellar fiux and background time-series 
with the full width at half maximum (FWHM) of the 
point response function (Demory et al. 2012a). We find 
for both AORs the lowest BIC to correspond to a model 
including a second-order position-dependent polynomial. 



4 



a second-order logarithmic ramp, and a time-dependent 
linear trend. Our analysis yields an RMS of 362 ppm and 
369 ppm per 5 min interval in the first and second AOR, 
respectively, with negligible contribution from correlated 
noise. 

3.1.2. Determination of the Stellar Density 

We perform a combined MCMC fit including our two 
Spzteer transits and the 61 HARPS radial velocities (RV) 
published in the discovery paper (Bonfils et al. 2012). 
The main goal of this step is to derive the stellar density 
from the Spitzer photometry (Seager & Mallen-Ornelas 
2003), to enable the derivation of the stellar and plane- 
tary physical parameters. The following system param- 
eters ('jump parameters') are left free in the MCMC 
fit, using uniform priors: the orbital period P, transit 
depth dF (planet-to-star area ratio, (_Rp/i?*)^). transit 
duration W , time of minimum light Tq, impact parame- 
ter b ~ a cos i/R^, the parameter K' = — e^P^/^, 
where K is the radial- velocity semi-amplitude, y^cosw 
and y^sincj. We use a quadratic law for the limb- 
darkening. We draw the theoretical values and corre- 
sponding uncertainties of the coefhcients ui and U2 from 
the tables of Claret & Bloemen (2011) for the TeS, \ogg, 
and [Fe/H] determinations reported in Section 4. We 
use the resulting distributions for ui and U2 as normal 
priors in our MCMC fit. We use the linear combina- 
tions ci = 2ui + U2 and C2 = mi — 2?i2 as jump param- 
eters, rather than ui and U2, to minimize the correla- 
tions of the resulting uncertainties (Holman et al. 2006). 
At each step of the MCMC fit, the stellar density is 
derived from this set of parameters and Kepler's third 
law. We run two chains of 10^ steps each, where the 
first 20% are discarded. We assess the good convergence 
and mixing of the chains employing the Gelman-Rubin 
statistic (Gelman & Rubin 1992). We add a 2.0 ms'^ 
jitter contribution in quadrature to the radial-velocity 
error bars to match the RMS of the residuals. This 
first combined run yields an eccentricity signal compat- 
ible with a circular orbit (-yecosw = —0.09 ± 0.14 and 
-Jesmoj = 0.00 ± 0.22). We therefore repeat the fit set- 
ting ^Je cos uj and ^/esmiij to zero. The difference in BIC 
between an eccentric and a circular orbit is ABIC = 15, 
translating to an odds ratio of ^-^1800, hence favoring 
the circular model we adopt in the following. Our anal- 
ysis assuming a circular orbit yields a stellar density 
Pi, = 2.911q 33P0 that we use as a constraint for the 
derivation of the stellar parameters in Sect. 4. The phase- 
folded Spitzer light curve is shown in Figure 3. 

3.2. Spectroscopic measurements 
3.2.1. FIRE Spectral Analysis 

We measure a number of spectral features in our FIRE 
spectrum of G J 3470 for the purpose of deriving its spec- 
troscopic properties, particularly [Fe/H] and Toff, us- 
ing recent calibrations presented for M dwarfs. We fol- 
low closely the prescriptions of Rojas-Ayala et al. (2012), 
Terrien et al. (2012), and Mann et al. (2012) for measur- 
ing equivalent widths (EWs) as well as water indices, and 
to estimate the pseudo-continuum needed for the above 
metallicity calibrations. 

For the Rojas-Ayala et al. (2012) calibration, we mea- 
sure the Na doublet (2.206 and 2.209 ^m) and the Ca 
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Fig . 3 . — G J 3470 b detrended and phase- folded transit light 
curve combining our two 4.5 fira Spitzer/TRAC visits, with the 
best-fit transit model superimposed (see Sect. 5). Data points 
arc binned in 2-minute intervals, and residuals are shown in the 
bottom panel. 



TABLE 1 

Metallicity Estimates [dex] for GJ 3470 from 
Near-Infrared and Visible Spectroscopy. 



Calibration Reference 


Bandpass 


Value 


Rojas-Ayala et al. (2012) 


K 


-1-0.15 ±0.17 


Terrien et al. (2012) 


H 


-1-0.25 ±0.12 


Terrien et al. (2012) 


K 


+0.19 ±0.12 


Mann et al. (2012) 


K 


+0.32 ±0.11 


Neves et al. (2012b) 


Visible 


+0.08 + 0.10 


Adopted 




+0.20 + 0.10 



triplet (2.261, 2.263, and 2.265 /im) following the integra- 
tion limits and continuum points of their Table 2. The 
pseudo-continuum flux of each feature is taken from a 
linear fit to the median flux within a 3 nm region around 
each continuum point. The water index, H2O-K2, is mea- 
sured following Eq. 5 of Rojas-Ayala et al. (2012). For 
the application of the H-hand and if -band metallicity 
calibrations of Terrien et al. (2012), we measure the EW 
of the Na (2.2074 ^m), Ca (1.6159, 1.6203, 2.2638 /.tm), 
and K (1.5171 pjn) features following the prescription de- 
tailed in Sect. 3.1 of their study. The pseudo-continuum 
is estimated by fitting a fourth order Legendre polyno- 
mial to the regions shown in their Fig. 1(A) for the H 
band, and in their Fig. 1(B) for the K band. The water 
indices, H2O-H and H2O-K, are also measured following 
the definitions in their paper. For the Mann et al. (2012) 
calibration, the EW of the metal-sensitive features F19 
(2.2079 /im), F20 (2.3242 /xm), and F22 (2.3844 /xm) in 
the A'-band are measured using the parameters listed in 
their Table 5. The pseudo-continuum is measured by a 
linear fit in the spectral regions specified in their Table 4, 
immediately redward and blueward of each feature. The 
water index used for this AT-band calibration is the same 
as the one described by Rojas-Ayala et al. (2012). 

The equivalent widths and spectral indices computed 
from our FIRE spectrum as described above yield the 
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metallicities for G J 3470 shown in the first four hues of 
Table 1. A spectroscopic estimate of the effective tem- 
perature of GJ3470 is obtained using the temperature- 
sensitive H2O-K2 index in the K band as defined by 
Rojas-Ayala et aL (2012). The result is 3750 ± 300 K. 

3.2.2. HARPS Spectral Analysis 

An additional spectroscopic estimate of the metallicity 
of GJ3470 is obtained from a recent calibration (Neves 
et al., in prep.) based on the visible-light HARPS spec- 
tra of Bonfils et al. (2012). This calibration was estab- 
lished on the basis of equivalent widths measured for a 
total of 4441 lines in the spectra of 55 stars from the 
HARPS volume-limited M-dwarf sample (Bonfils et al. 
2013), and is anchored on existing photometric cali- 
brations for metallicity (Neves et al. 2012a) and effec- 
tive temperature (Casagrande et al. 2008). The proce- 
dure, described briefiy in the Appendix of the study by 
Neves et al. (2012a), achieves an improved precision over 
previous methods of 0.10 dex. The result of this measure- 
ment for GJ 3470 is +0.08 ± 0.10, and is collected also in 
Table 1. 

The five estimates of [Fe/H] from the FIRE and 
HARPS spectra are consistent with each other, and 
we therefore adopt for the remainder of the paper the 
weighted average, [Fe/H] = +0.20 ± 0.10, in which the 
uncertainty is a more conservative estimate than the for- 
mal error of the mean. 

4. STELLAR CHARACTERIZATION 

Mass (M^) and radius (i?j,) estimates for exoplanet 
host stars are typically obtained by appealing to stel- 
lar evolution models. For M dwarfs this, too, has 
been problematic (beyond the challenges for determin- 
ing Toff and [Fe/H] alluded to earlier in Sect. 2.3) be- 
cause of known disagreements between predictions from 
theory and accurate measurements of M^, and i?^ for 
low-mass stars in double-lined eclipsing binaries (see, 
e.g., Torres 2012, and references therein). We there- 
fore rely here exclusively on empirical relations, on the 
mean stellar density inferred from our Spitzer light curve 
in Sect. 3.1.2 (p* = 2.9ltQ;33 ^0)1 and on brightness 
measurements for GJ 3470 from 2MASS and in the op- 
tical (V = 12.33 ± 0.01; Weis 1986; Evans et al. 2002; 
Zacharias et al. 2013). For a given parallax and ignoring 
extinction, the near-infrared mass-luminosity (M-L) re- 
lations of Delfosse et al. (2000) provide estimates of the 
absolute mass, and are insensitive to metallicity. On 
the other hand, the surface-brightness (SB) relations by 
Kervella et al. (2004) allow one to estimate the angu- 
lar diameter, which may be converted to a linear ra- 
dius with knowledge of the parallax. The latter rela- 
tions are valid for [Fe/H] between —0.5 and +0.5, and 
are thus applicable to G J 3470, with its metallicity of 
[Fe/H] = +0.20 ± 0.10. While a trigonometric paral- 
lax has not been measured for this star, we may use the 
above relations simultaneously to solve for the distance 
that yields values of M^, and i?^ consistent with the mea- 
sured mean density. 

We proceeded in a Monte Carlo fashion, drawing all 
measured quantities {VJHKg photometry and Spitzer 
mean density) from appropriate Gaussian distributions. 
For each set of draws we solve for the value of the par- 
allax that gives a mass and radius through the M-L 



and SB relations resulting in a mean density equal to 
the randomly drawn value of pi, for the set. We re- 
peat the process 10^ times, and adopt as final values 
the mode of the corresponding posterior probability dis- 
tributions, assigning l-cr uncertainties given by the 15.85 
and 84.13 percentiles of those distributions. We obtain 
A'U = 0.539lEJ |^4^ M© and = 0.568tJJ:[J^^ i?©, and a 
parallax of tt = 32.4lt^ g mas, corresponding to a distance 
of 30.7^17 pc. The mass is an average of the J-, 
and if -band relations by Delfosse et al. (2000), each of 
which is assumed conservatively to carry an uncertainty 
of 10%. The radius is an average of the two surface- 
brightness relations of Kervella et al. (2004) that yield 
the smallest scatter in the angular diameters (about 1% 
for the relations that depend on V~H and V~K). Prior 
to using them, the 2MASS magnitudes are converted to 
the native photometric system of the M-L and SB rela- 
tions (CIT and Johnson, respectively) using the trans- 
formations of Carpenter (2001). The uncertainties listed 
above include all photometric errors, the error in p^,, as 
well as the scatter of the empirical relations. We note 
that our stellar mass is very close to that reported by 
Bonfils et al. (2012), but our radius is 13% larger. 

As a check on the above absolute radius determination, 
we obtain additional estimates of i?* from color indices 
and the calibrations recently published by Boyajian et al. 
(2012), which are based on angular diameter measure- 
ments from the CHARA interferometer and HIPPAR- 
COS parallaxes, and have a dependence on metallicity. 
Results using V — J ^ V~H, and V — Kg for the mea- 
sured metallicity of GJ 3470 give very consistent values 
for Ri, averaging 0.513 + 0.043 i?©, in agreement with our 
Spitzer-hased determination within about la. Figure 4 
displays the location of GJ 3470 in the mass-radius dia- 
gram for low-mass stars, along with the measurements 
for all other such objects in double-lined eclipsing bina- 
ries that have relative measurement precisions under 5% 
for Af^ and i?*. The constraint afforded by the mean 
stellar density is also indicated. 

While an estimate of the effective temperature of 
the star was obtained earlier using our FIRE spec- 
trum, the precision is relatively low. We obtain a fur- 
ther estimate using the color/temperature calibrations 
of Boyajian et al. (2012), which are based on bolometric 
fluxes and angular diameter measurements, and include 
metallicity terms. The V — J, V — H, and V — Kg in- 
dices along with our adopted value of [Fe/H] lead to a 
weighted average temperature of 3630 ± 100 K. A final 
Tcff estimate is inferred from the same three indices and 
the color/temperature calibrations of Casagrande et al. 
(2008), which rely on the Infrared Flux Method. How- 
ever, these relations do not take into account the metal- 
licity, and implicitly assume a composition near solar 
whereas G J 3470 is metal-rich. Therefore, the resulting 
estimate (3360+ 100 K) requires an adjustment for metal- 
licity. We determine this by using the stellar evolution 
models of Dotter et al. (2008) in a differential sense, first 
reading off from a [Fe/H] = +0.20 isochrone the stellar 
mass that yields the same color indices as we measure, 
and then comparing the corresponding temperature with 
that for a star of the same mass on a solar-metallicity 
isochrone. This exercise is insensitive to the age adopted 
for the isochrone. Consistent results using the three color 
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Fig. 4. — Mass and radius estimates for GJ 3470 (large red dot) 
compared with measurements for other low-mass stars in double- 
lined eclipsing binaries with relative errors in M^, and R^, loss than 
5% (see Torres 2012). The constraint on the mean stellar density 
derived from our Spitzer observations is indicated with the dashed 
line (dotted lines representing the 1-cr uncertainties). 

indices separately give an average correction of -1-140 K, 
which resuhs in a final temperature of 3500 ± 150 K. As 
the two photometric determinations above are consistent 
with each other and with the spectroscopic determination 
in Sect. 3.2.1, we adopt the weighted average of the three 
values, Teff = 3600 ± 100 K. 

While this paper was under review, we learned that 
Pineda et al. (2013) performed an independent charac- 
terization of GJ3470's stellar properties. We refer the 
reader to that study for a description of their analysis 
and results. 

5. PLANETARY AND ORBITAL PARAMETERS 

Our final parameters for G J 3470 b were derived using 
the stellar properties from the preceding section, and a 
MCMC fit analogous to that described in Sect. 3 with 
the addition of a prior on the stellar mass. For this prior 
we used the posterior probability distribution derived in 
Sect. 4, drawing from it a random value of the mass 
at each step of the MCMC fit. As before, we included 
the light curves from both Spitzer visits, along with the 
61 HARPS radial velocities reported by Bonfils ct al. 
(2012). The results are presented in Table 2, where the 
value adopted for each parameter is the median of the 
corresponding marginalized posterior distribution from 
the MCMC fit. Error bars are the corresponding 68.3% 
probability intervals from the same distributions. The 
final model and phase-folded Spitzer light curves are dis- 
played in Figure 3. 

We find for GJ 3470b a radius of Rp = 4.83t°J2 R^, 
which is 13% larger than previously reported in the lit- 
erature. This increase is driven mainly by the larger 
stellar radius from Sect. 4. Combining the planetary ra- 
dius with the mass Alp — 13.9't\'\MQ that relics on the 
radial- velocity data set yields a very low planetary den- 
sity of pp — 0.72^0 12 gcm~^, which is 33% smaller than 
the estimate in the discovery paper. These planetary 
parameters arc also listed in Table 2. 

Finally, we performed a new fit for the purpose of as- 



sessing the robustness of the orbital period determination 
for GJ 3470 b, which in our solution is constrained both 
by the two 5piteer transits and the radial velocities. How- 
ever, the two Spitzer Yisits are consecutive (2012 June 11 
and 15), so the lever arm for the orbital period determi- 
nation is very short. We therefore incorporated the two 
TRAPPIST transit light curves from Bonfils et al. (2012), 
as well as the ones from EulerCam and the NITES tele- 
scope. The light curves from the first two sources show 
only the ingress portion of the transit, but may still be 
combined with our two full Spitzer light curves that con- 
strain the transit shape, if we assume the latter does 
not change across wavelengths. The NITES light curve 
has a higher level of correlated noise, but does cover the 
transit completely. These ground-based light curves were 
obtained between February and April of 2012, and there- 
fore contribute to build up a much longer baseline. 

As expected, most of the system parameters in this 
new fit are tightly constrained by the Spitzer photometry 
alone, but the period is considerably improved. The new 
value is included in Table 2, and is only 19 ± 11 s shorter 
than the one that relics on the two 5piteer transits alone. 

6. INTERIOR COMPOSITION OF A LOW-DENSITY 
EXO-NEPTUNE 

G J 3470 b presents a valuable test case for planet for- 
mation and evolution theories. It stands out from the 
crowd of accumulating transiting exo-Ncptunes due to 
its low mean density and bright M dwarf host star. 
GJ3470b's measured radius is 20% ± 6% larger than 
Uranus (i?g = 4.01 R^) despite having a similar mass 
(Mg = 14.5 Mq). The planet radius corresponds to 
roughly 20% of its Roche-lobe radius. Among currently 
known low-mass [Mp < 30 M^) transiting planets, only 
the Kepler-11, Kepler-18, Kepler-30, and HAT-P-26 sys- 
tems have planets with lower densities (see Figure 5). 

G J 3470 b must have acquired H/Hc gas from the pro- 
toplanetary nebula. Alternative gas layer sources such as 
sublimated ices and outgassing from a rocky interior may 
be important for less massive, more dense planets (such 
as GJ 1214b; Rogers & Seagcr 2010b), but cannot be the 
dominant gas layer source for GJ 3470 b. Its bulk density 
is too low for astrophysical ices (H2O, CO2, etc.) alone 
to comprise the planet volatiles; significant quantities of 
light gases (hydrogen and helium) must be present. Fur- 
ther, GJ 3470 b's gas layer is too voluminous to have been 
formed by outgassing of light gases during formation; the 
planet radius exceeds the upper limit for outgassed plan- 
ets from Rogers ct al. (2011). 

Nebular H/He contributes between 5% and 24% to 
GJ 3470 b's mass, according to our interior structure 
models. Following Rogers & Seagcr (2010b), we apply a 
fully differentiated model for the planet's interior struc- 
ture consisting of (from the center of the planet outward) 
an iron core, silicate layer, ice layer, and H/He gas en- 
velope to explore which bulk compositions are consistent 
with the measured mass and radius of G J 3470 b. Both 
the planet's bond albedo A (which scales the equilibrium 

temperature Toq = (1 - A)^/^ (683 ± 27) K), and the 
planet's intrinsic luminosity Lp (a proxy for the poorly 
constrained age of the planet) are unknown. We adopt 
fiducial values of ^ = 0.3 and Lp/Mp = 10"!° Wkg-\ 
while also exploring the ranges of A = to 0.5 and 
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TABLE 2 

Adopted system parameters for GJ 3470 from our 
MCMC FIT OF Sect. 5. 



Parameter 



Value 



Jump parameters 



Planet/star area ratio Rp/Rs 

b = acosi/Rf, [R^,] 

Transit width W [d] 

To - 2,450,000 [BJDtdb] 

Orbital period P [d]^ 

RV K' [ms-idi/^*] 

y/e cos ui 

ci = 2ui + U2 
C2 = ui — 2u2 

Stellar parameters 



O.U7798_o 00045 

0.0791 ± 0.0005 
6090.47701 ± 0.00010 
3.33665 ± 0.00005 
13.4 ± 1.2 
0.0 (fixed) 
0.0 (fixed) 
0.246 ± 0.027 
-0.329 ± 0.020 



Ml 
"2 

Mean density p* [pq] 
Surface gravity log [cgs] 
Mass M* [Mq]^ 
Radius i?* [Rq\^ 
Parallax tt [mas]'' 
Distance [pc]'' 

Effective temperature T^g [K]'' 
Metallicity [Fe/H] [dex]'' 

Planetary parameters 



0.033 ±0.015 
0.181 ± 0.010 

4.658 ± 0.035 

U.003_o 

U.ODf5_o 031 
1-2.1 



1.9 
2.1 



32.4] 
30.7"_i 7 
3600 ± 100 
+0.20 ± 0.10 



RV semi-amplitude K \ms ^] 
Orbital semi-major axis a [AU] 
Orbital inclination i [deg] 
Mean density pp [g cm"^] 
Surface gravity log [cgs] 
Mass Mp [M®] 
Radius Rp [R^\ 

Individual transit timings 



8.9 ± 1.1 



0.03557; 
88.3 
0.72 
2.76 
13.9 
4.83 



+0.00096 
-0.00100 
+0.5 
0.4 
+0.13 
0.12 
+0.06 
0.07 
-1.5 
-1.4 
+0.22 
0.21 



To 1 - 2,450,000 
To,2 - 2,450,000 



BJDtdb: 
BJDtdb. 



6090.47705 ± 0.00014 
6093.81372 ± 0.00015 



^ Derived using our two Spitzer light curves along with pub- 
hshed ground based photometry and RVs (see Sect. 5). 
^ Parameters derived either in Sect. 3 or in Sect. 4, and 
repeated here for convenience. 



Lp/Mp = Wkg-i to 10-9-5 Wkg"i. Figure 6 

presents the H/Hc gas mass fraction (Mxy/Mp) in our 
models as a function of the Fe-silieate-H20 abundances of 
the heavy element interior (assuming the median values 
of the planet mass and radius, and our nominal planet 
energy budget parameters). Varying the planet mass and 
radius within their l-tr bounds, and considering a range 
of plausible planet energy budgets affects the H/He mass 
fractions by up to ±0.05. For a rocky Earth-like heavy 
element interior composition (32% Fe, 68% silicate, 0% 
H2O), GJ3470b's H/He envelope mass is constrained to 
Mxv/Mp = 0.16±0.05, while for a denser iron-enhanced 
Mercury-like rocky interior (70% Fe, 30% silicate, 0% 
H2O), Mxv/Mp = 0.17 ± 0.05. Less H/He is needed 
if G J 3470 b has an ice- rich interior composition; for in- 
stance, for a heavy element interior with 16% Fe, 34% 
silicate, 50% H2O, Mxv/Mp = 0.12l° °|. 

Which heavy element interior compositions arc plausi- 
ble for G J 3470 b? The planet interior icc-to-rock ratio 
is not constrained by measurements of the planets mass 




15 20 

M {M ) 

Fig. 5. — Mass-radius relationships of small transiting planets. 
G J 3470 b is highlighted in red. Other small transiting oxoplanets 
with dynamical mass constraints (CoRoT-7b, Kepler-4b, Kepler- 
10 b, Kepler-llb, c, d, e, f, Kepler-18b, c, d, Kepler-19b, Kepler- 
20 b, c, d, Kepler-30b, d, Kepler-36b, c, 55Cnce, GJ 1214 b, 
GJ436b, HAT-P-llb, and HAT-P-26b) are plotted in black. The 
solar system planets are indicated with solid triangles. The curves 
are illustrative constant-temperature mass-radius relations from 
Seager et al. (2007). The solid lines are homogeneous-composition 
planets: water ice (blue solid), MgSiOs perovskite (red solid), and 
iron (magenta solid). The non-solid lines are mass-radius relations 
for differentiated planets: 75% water ice, 22% silicate shell, and 3% 
iron core (blue dashed); Ganymede-like with 45% water ice, 48.5% 
silicate shell, and 6.5% iron core (blue dot-dashed); 25% water ice, 
52.5% silicate shell, and 22.5% iron core (blue dotted); Earth-like 
with 67.5% silicate mantle and 32.5% iron core (red dashed); and 
Mercury-like with 30% silicate mantle and 70% iron core. 

and radius alone, so we look to planet formation the- 
ory for insights. If G J 3470 b formed beyond the snow 
line and migrated inwards to its current orbit, its heavy 
element interior would be ice-rich. If instead G J 3470 b 
formed in situ (inside the snow line) its heavy element 
interior would be rock dominated with a lower propor- 
tion of ices. Theoretical predictions for how much ice is 
likely included in planets formed inside the snow line of 
M dwarfs are a topic of ongoing debate. Ogihara & Ida 
(2008) proposed that migration of planetesimals from be- 
yond the snow line could supply icy material to the in- 
ner regions of the protoplanetary disk. On the other 
hand, Lissauer (2007) and Kennedy et al. (2007) predict 
that planets and planetesimals formed within 1 AU of 
M dwarfs are unlikely to have large volatile inventories 
when the effect of the M dwarfs' pre-main sequence lu- 
minosity evolution is taken into account. In Figure 6 we 
present interior bulk compositions for the full range of 
ice-to-rock ratios. 

7. SUMMARY 

Our 4.5 /im Spitzer observations have enabled us to re- 
fine the planetary and system parameters of the Neptune- 
size planet G J 3470 b, improving its radius to Rp = 
4.8±O.2i?0, which is 13% larger than previously reported 
in the literature. As a result, the revised planetary den- 
sity, Pp = 0.72±0.13 gcm"^, is 33% smaller than before. 
These changes come mostly from revisions of the stellar 
parameters (particularly R^), which have been frustrat- 
ingly difficult to determine accurately in the past due to 
known discrepancies between observations and standard 
stellar evolution models for lower main-sequence stars. 
In this paper we have relied for this only on empiri- 
cal mass-luminosity and surface brightness relations that 
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Fig. 6. — Fraction of GJ 3470 b's mass contributed by H/He, as a 
function of the planet's heavy-element interior composition. Each 
point within the diagram corresponds to a specific combination of 
Fe, Mgo.gFeo.iSiOs, and H2O (by mass) in the heavy element inte- 
rior of GJ 3470 b. For instructions on how to read ternary diagrams 
see, e.g., Valencia ct al. (2007); Zcng & Seager (2008). Each con- 
tour is labeled with Mxy/^^p for ouv fiducial model parameters 
(median Mp, median Rp, A = 0.3, and Lp/Mp = 10'^° Wkg-^). 
Uncertainties in the planet mass, radius, and energy budget can 
affect Mxy/Mp by ±0.03 to 0.05. 



have been widely employed in other contexts, and on the 
strong constraint on the mean stellar density provided by 
our Spitzer observations. In the process we have inferred 
an accurate distance for the star. 

G J 3470 b provides a valuable example of an extremely 
low-density planet, representative of a significant portion 
of the exoplanet candidates found by the Kepler mission 
to date. The brightness of the host star {Kg = 7.99) com- 
bined with its large planet-to-star radius ratio renders 
G J 3470 b a promising candidate for future atmospheric 
characterization, which could provide clues on its forma- 
tion pathway. Indeed, GJ 3470 b's low surface gravity 
translates to a large atmospheric scale height for a given 
atmospheric composition, favoring follow-up studies ap- 
plying transmission spectroscopy. GJ 3470 b, G J 436 b, 
and G J 1214 b are a remarkable sample of volatile-rich 



planets orbiting bright nearby stars, pushing the field 
of comparative exoplanetology further towards low-mass 
planets. 
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